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ABSTRACT 

: 

It has been shown, in the case of meson photoproduction, that the power-law falloff of these 
reactions can be described by lowest order (real) sum rules, at moderate momentum transfer. 
The phases of these processes, in this regime, are usually thought to be non-perturbative. In 
a sum rule framework, however, they can possibly be described by radiative corrections to the 
^ . hadronic spectral densities of the corresponding helicities, which become complex functions 
'. to order a^, and the effects of interference can be strongly enhanced by the presence of the 
^ I vacuum condensates in the dispersion relation. It is shown that the imaginary parts of these 
if^ ' complex corrections have a factorized form and can be evaluated in a systematic fashion, while 
O ■ their real parts, at the same perturbative order, are down by at least 2 powers of momentum 
Q>^ . transfer. The analysis is done at two loop level, combining dimensional regularization and 
"Th ! light-cone methods. The calculations are performed for all the independent set of scalar 
PhI diagrams generated by the OPE. The analytical bounds are identified and discussed. 
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1 Introduction 



Sum rules have been a valuable tool in the investigation of the intermediate energy region 
of QCD. More recently, it has been suggested that such methods can find application in the 
context of elastic scatterings of Compton type 0. Compared to sum rules for form factors, 
which involve dispersion relations only at fixed t, the description of elastic scattering has 
been performed by sum rule based on dispersion relations at both s and t fixed, and with a 
limited domain of analyticity in the remaining variables. For scatterings at intermediate angles 
such new description of 2-photon elastic processes compares favourably well with standard 
factorization theorems 0]. Such extended sum rule methods can be helpful in the attempt to 
understand on a firmer ground the transition to perturbative QCD in simple elastic processes, 
and can have potential application to a new entire class of direct reactions of Compton type 
and to the corresponding crossed channels. 

In fact, two-photon collisions are among the most interesting ways to study these important 
aspects of elastic scattering at moderate momentum and energy transfers. 

Together with Li, we have analized in detail the main features of the sum rules for pion 
Compton scattering and pion photoproduction [|16| and compared their predictions with those 
derived from modified and standard factorization theorems fl^ . 



All the work presented so far |^ for these reactions is based on real sum rules for the 
two helicities Hi and H2. 

The analysis of the radiative corrections to the lowest order result are, however, of crucial 
importance in the study of the phases of these amplitudes, since nontrivial interference effects 
between those can be generated by such corrections. Differently from the usual prescriptions 
based on factorization theorems, such effects, in the sum rule context, can be strongly enhanced 
due to the presence of the vacuum condensates as fundamental parameters in the dispersion 
relations. 

It is widely believed that in the intermediate energy region of QCD (Q^ ~ 4-7 Gev"^) such 
phases cannot be incorporated in a direct perturbative treatment and are of non-perturbative 
origin. In the sum rule approach interference effects are generated by complex contributions 
(due to gluonic exchange) to the lowest order (real) spectral densities. 

The evaluation of such corrections, however, is a notrivial task and requires considerable 
effort. In this work (section 3) we discuss in great detail the properties of analyticity of the 
lowest order spectral densities of those correlators which interpolate with processes of Compton 
type 0. The motivation of our study relies mainly on the fact that the region of analyticity 
of 4-point correlators, compared to that for the 3-point correlators which are commonly used 
in the sum rules for vertex functions, has unobvious features. Such a region turns out to 
be bound in size by w-channel singularities, a moving w-channel cut, whose location varies 
when the other variables which appear in the dispersion integral are varied. The study of the 
properties of analyticity of the spectral densities which appear in sum rules of these type is 
closely related to the mass parameters in the correlator of the 4-currents from which the OPE 
is calculated and the sum rules are derived. 



2 



A preliminary discussion of the mass dependence has been presented in ref. 0, where 
power corrections to the sum rule for a specific combination of the helicities of pion Compton 
scattering have been presented. An artificial mass dependence is generated in the calculation 
of such corrections (this point is briefly illustrated in section 4) and we have argued, in the 
same work, that this feature was not going to affect the size of the analyticity region or the 
validity of the sum rule, for small values of the mass. 

Here we are going to give a general proof of these statements together with an explicit 
description of the analytical bound for the dispersion relation, generated by a finite mass 
dependence in the expansion of the correlator. 

Such thresholds, absent in the case of form factors, cannot be ignored in the case of s and t 
dependent dispersion relations. This bound is given by a quadratic equation, and shows up as 
a singularity of the lowest order spectral density (the vanishing of its denominatror). We show 
in sections 3 and 4 that moderate values of s and t keep the bound far away from the (finite) 
region in which the OPE is calculated. A similar (mass-dependent) bound appears for the 
diagrams describing the power corrections. The presence of these bounds is justified by the 
fact that each coefficient of the OPE has leading singularities (also called Landau surfaces) 
generated when all the internal lines in the perturbative expansion go on shell. For Compton 
scattering these surfaces corresponds to a m channel threshold. We show, however, that these 
singularities, though important in a more general analysis (in mass dependent correlators), 
disappear in the massless limit. Therefore, in this particular case, a simplified structure of the 
OPE emerges, which is fully exploited in the calculation of the radiative corrections. We are 
also able to write down a dispersive representation of the main integrals (the coefficients of 
the OPE) which is unbound in the plane of the two dispersive variables. The proof is limited 
to the lowest order result. 

In this simplified case we prove that a limited applicability of Borel methods [||, H (Ap- 
pendix C) is possible, and that the leading spectral density can be re-obtained elegantly 
without using the usual Cutkosky rules. This is consistent with our results of section 3. 

The absence of these w-channel cuts (in the massless case) in the lowest order correlator 
is crucial for the evaluation of the radiative corrections. In fact we are looking for imaginary 
parts to the spectral densities and, to next order in perturbation theory, these cuts reappear 
in specific subdiagrams. The analysis presented in sections 3 and 4 is therefore preliminary 
to our subsequent discussion. In section 5 we illustrate, using various results of dispersion 
theory, how to organize the calculation of the radiative corrections to the spectral densities 
in a systematic way. The analysis of the complex parts in the diagrammatic expansion of the 
spectral density, is illustrated for all the independent set of diagrams, and a simple factorized 
structure (in terms of specific sub — cuts) of the OPE emerges. 

Results of dispersion theory to lower order are then used to decide of the real or complex 
behaviour of such functions. Remarkably, thanks to this simple factorized feature, we show 
that such corrections can be evaluatedd in closed form. For this purpose, we combine dimen- 
sional regularization and light-cone methods, and show that the divergences in the OPE can 
be isolated as poles in e = n — 4. 
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Our conclusions are in section 6. Appendices A-C cover technical derivations of sections 
3-5, while Appendix D illustrates, in a self contained way, how to extend Borel methods 
to Compton scattering. 



2 Spectral densities for 4-point correlators 

Sum Rules relate the timelike region of a correlator to its spacelike part by a dispersion relation 
^, |12|. In the case of pion Compton scattering, for instance, the 4-point correlator with 
non vanishing projection over the invariant amplitudes (the helicities) of the process is given 
by i 

^a^,uxipl,pl, s,t) = i J d'^xd^yd^zexp{-ipi ■ x + ip2 ■ y - iqi ■ z) 

x{0\T{r]^{y)J,{z)M0)4{x))\0) , (1) 

where 

2_ 1 - 

J^l = ^U'-i^u - -d'-i^d, 7]^ = u-f5'jad (2) 

are the electromagnetic and axial currents, respectively, of up and down quarks, qi and q2 are 
on shell moment carried by the two physically polarized photons (see Fig. 1). 

The two pion momenta are denoted by pi and p2, with si = pi and S2 = pi being the 
virtualities. 

The invariant amplitudes of pion Compton scattering are then extracted from the physical 
expansion of this correlator for timelike momenta pi,P2- The matrix element which interpo- 
lates with pion states is given by 

M,x = ^| d'ye~''^^y{p2\T{My)J,{0)) \p,) , (3) 

and can be related to the two helicities Hi , H2 as 

M^^ = ifi(s, t)e(^)^e(^)^ + H2is, t)e^^^^e^^^'', (4) 

where e'-^-' and e*^^-* are helicity vectors defined in @, §]. 

The spectral densities appearing in the sum rules are characterized by a resonance con- 
tribution at the pion pole A'''^'*(si, S2, s, t), and by a continuum one Ap'^''*(si, S2, s, t), at large 
values of the virtualities si and S2 0, 0]- A dispersion relation than allows us to connect the 
spacelike region of the 4-current correlator in eq. (0) to its timelike behaviour in p1,P2 at fixed 
angle —t/s. 

The resulting sum rule for Hi{s,t) {i = 1, 2) can be written down in terms of a local 
duality contribution and of power corrections of quarks and gluons as 
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t 

/ dSi dS2P^''' + —{G^) dSi dS2 

Jo Jo vr JO Jo 

+Cr'''7ra.mf) . (5) 

where sq, the local duality interval, characterizes the resonant region of the correlator. 
The spectral densities pP^rt^ pquark^ ^giuon calculated by the OPE. They are polynomial in 
a variable Q{si, S2,t) which has a natural interpretation in the light-cone frame of the pion' 
lines (see Appendix A). 

The evaluation of all the lowest dimesnional contributions to the OPE of the two helicities 



has been presented in ref. |16 



For instance, the covariant expression of pf is of the form 

pert _ _5 RjjQ^, S, Si, S2) . X 

~ 24772 g4(2Q2 _ 3)2(4Q4 _ SiS2yi2Q^S - SiS2y ' ^ ' 



where 
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RiiQ^, S, Si, S2) = Si, S2)Q^". (7) 

n=0 



The explicit forms of the coefficients ai^n can be found in [|T^]. The contributions from the 
lowest dimensional vacuum condensates are also given in that work. 

As shown below in eq. (^), the expression of Q{si,S2,t) is non trivial and although it 
can be viewed as the basic momentum scale which characterizes the power law fall off of 
the spectral densities (see 0), (which are simple rational functions when expressed in this 
variable), it is however not exactly on the same ground as s and t in characterizing the 
"dimensional counting" of the sum rules (as it is usually done in standard factorization [P^). 
There are two main reasons for this: 1) Q contains the sum of both s and |t| at larger energies 
(which in our case are supposed to be of the order of 5-7 Gev'^ or so @]) 2) it appears as a 
variable (under the local duality interval) and therefore it is not just a fixed external invariant. 



The power law falloff of the spectral density, however, shows that [y, |T^, generically 



p{Si,S2,S,t) ^^ + ^,+... (8) 

where Ci , C2, etc, are complicated functions of the pion virtualities and contain the input from 
the condensates. Therefore, contributions from "higher-twists" are embodied in an obvious 
form by the dispersion relation. There are some specifc features of this approach which clearly 
differentiate it from the perturbative picture: 1) to lowest order, no gluons are involved, 2) the 
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power law falloff is much more suppressed (a coherence effect) compared to the perturbative 
one and 3) the sum rule prediction starts dominating compared to the perturbative one toward 
the region of forward scattering. Features 1) and 2) are known from the literature on sum 
rules for form factors. However, in this latter case, the contributions from the condensates 
grows as a power of " Q" and soon overtakes the (decaying) contribution from the perturbative 
coefficient of the OPE. In our case, instead, a suppression is generated -for scattering at fixed 
angle- at the photon vertex by a virtual quark line connecting the two photons and eliminates 
this growth. In fact , compared to form factors, a dependence on both s and t is now availiable 
in the sum rule. The shortcome of this, however, is that the stability analysis of the sum rules 



(in both s and t) becomes much more complex |jT6 



In order to describe, for instance, interference effects in photoproduction processes (7 7 — > 
7r+7r~) and in other related processes by a sum rule to order as, one needs the imaginary parts 
of the radiative corrections to the spectral densities. The real parts are also, as a matter of 
fact, needed, however an actual calculation shows that these are suppressed by a factor l/Q^ 
compared to the imaginary ones. The argument is quite simple and can be better understood 
from a reading of section 5. Basically, imaginary parts are generated by pinched diagrams 
containing 5-particle cuts, while real parts require only 4-particle cuts. This implies that more 
propagators of large virtuality (1/Q^) appear in the integrals describing real discontinuities, 
compared to those contributing to the imaginary ones. One can also define suitable observable 
were the contributions from the imaginary parts are enhanced compared to the real ones. 

Here is a simple example of how this may be happen. 

If we include radiative corrections. Hi and H2 are determined by the complex spectral 
densities 



Pi(si, S2, s, t) = pf^ + Us (pf ^ + i pf^) + pf (9) 

where the suffix i runs from 1 to 2. p^^ is a short notation for the contributions coming from 
the power corrections. It is not difficult to choose specific observable which are more sensitive 
to these complex parts than to the real ones. For example, simple counting arguments are 
sufficient to show that the difference 



cr = (j(7 7 ^ tttt) - o-(7ij 7/j ^ 7r+ vr ) (10) 

is dominated by the lowest order spectral density and by the purely imaginary parts of 
their radiative corrections (po- ~ P^i^ P2 + P^? pi"^)- we are going to show in the next few 
sections, the imaginary parts of the sum rules of meson photoproduction can be classified 
according to a well defined factorized structure and are calculable in a close form. 

Notice that the sum rule studied so far have been calculated in the approximation of 
massless quarks in the expansion of the correlator. Although this is clearly apppropriate in 
the case of the pion, may not be appropriate for the investigation of heavier mesons. 

In fact, the method can be easily applied to all the other reactions obtained by simple 
crossing (in the s, t and u channels) of the Compton scattering amplitude, and in particular. 
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to processes of the form 7 + 7 A + A, where A can be, in principle, any meson. 

The evaluation of the sum rule, for these more general reactions, encounters a special 
difficulty related to the presence of a mass m in the diagrammatic expansion. As we are going 
to point out in the next few sections, this problem does not show up in the case of the form 
factor since the spectral density - at least to lowest order - is mass independent 0. Such 
issues, although of technical nature, deserve special attention in the sum rule context. Our 
discussion, here, is ground-breaking since there are no available methods in the calculation of 
the OPE of the double spectral densities beyond 2- and 3-point functions. In particular we 
point out that the use of the Borel transform in the isolation of such spectral functions 
must be handled with a certain care when the rank of the correlator increases (such as for 2 
photon processes). We refer to the discussion in Appendix C for a presentation of the Borel 
method and its partial generalization to the evaluation of the lowest order spectral functions 
for Compton scattering. 

In the case of the form factor, the radiative corrections, even in the massive case, can 
be obtained by the Borel method |10] which is, in fact, the easiest way to evaluate such 



contributions. As we are going to discuss next, in Compton scattering, these methods, 
which rely on simple properties of the Borel transform and on the existence of unbound 
dispersion relations for the corresponding Feynman diagrams, are of limited use. Although 
the particular kinematics chosen in ref. for the analysis of this reaction by sum rules (the 
requirement of treating on-shell physical photons) brings our treatment closer to the one in 
the form factor case, yet new analytical features emerge in the leading order spectral densities 
of these processes and in their power corrections Therefore, we proceed in sections 4 and 
5 to develope our own approach to the identification of the complex contributions to the sum 
rules. 

Our discussion, in the next sections, rehes only on the analysis of massless and massive 
scalar amplitudes. There is no need, in fact, to consider the full expression of the relevant 
diagrams for processes of these type (and this is true also in the form factor) since the basic 
difficulties in the analytic properties of the spectral densities are all contained in the scalar 
amplitudes. 



3 Massive correlators 

Here we start our discussion of the role of a mass dependence by considering 3-point functions. 
The general features of the sum rule method for the analysis of form factors can be found 
in refs. [0 and |^. In this case the region of analyticity (for any 3-current correlator 
which interpolates with a specific form factor) is unbound and the dispersion integral extends 
up to infinity. In fact the basic spectral density is obtained from the integral (the triangle 
singularity) JIT], H 

A3(si, S2, t)= I d^k5+{e - ml)6+{{p, - ky - ml)5+{ip, -k + q^f - ml), (11) 
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The propagators of the cut hues have been replaced, as usual, by "positive" delta functions 



^2 _ ^2 



27Ti6+{k' -m'), (12) 



as prescribed by Cutkosky rules. The relevant diagram is depicted in Fig. lb. It gives the 
leading spectral function, for form factors, to lowest order in as |^ . The dashed lines describe 
the usual cutting rules for the propagators. The spectral density for such a diagram takes the 
form (with pf = Si, P2 = S2) 

^3(^1, ^2, t) = \ '\ (13) 

2(A(si,S2,t)) 

where 

A(si,S2,t) = ((si + S2-t)'-4siS2))'^'. (14) 

is the Mandelstam function. The evaluation of ([TI|) can be easily carried out in the Breit 
frame of the two (pion) lines characterized by the virtualities pi and p^. It is also easy to 
realize that ( pTSD does not develope any additional singularity for positive si, S2, and negative 
t {t = {p2 — PiY fixed). 

Another important observation, in the case of form factors, is the independence of the 
discontinuity eq. ([TT|) from any mass in the propagator (see also |11]). Eq. (11) is in fact 



regular for any si and S2 over the entire complex plane of each of these variables. This allows 
us to write down a double dispersive representation of the massless triangle diagram in a 
rather straightforward way as simply as 



Tz{pl,pl,t) 



k'^{pi - kY{p2 - k) 



2 



00 



dsi [ ds2-, ^' ^ ON + subtr. 

Jo [Sl - Pi)[S2 - P2) 



(15) 



In eq. (0) we have omitted subtractions and single dispersive contributions to this integral. 
These two features 1) mass independence of the double spectral function and 2) its regularity 
over the entire si and S2 complex planes, are not preserved in Compton scattering, except for 
the massless case. Additional care is therefore required when we switch on a mass dependence 
in the basic diagrams, since new thresholds related to the w-cut in the dispersion variables 
automatically appears in the complex planes of the two virtualities si and S2- 

To be specific, let's consider the full contribution to the box diagram with a momentum 
flow chosen as in Fig. lb, and for simplicity, let's restrict our considerations to the scalar case. 
We consider the following 4-point function 

j ^^g^ 

J k'^ipi — ky(pi — k + qiY{p2 — k)"^ 
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with massless propagators. At fixed angle and with 

s + t + u = pl+pl s>0, t<0, n<0 (17) 

and moderately large s, t, u, assuming the existence of a region of analyticity of moderate size 
in si and S2 (see Fig. 2), we can write down a spectral representation for such integral of the 
form 



An^ Jji J^2 {si - pj){s2 - P2) 

where the contours 71^2, both of a radius A^, are again chosen as in Fig. 2. If we introduce a 
double spectral function A(si, S2, s, t), we can rewrite T4 as 

2 ,^ ^ A{s^,S2,S,t) 

T4{Pi,P2,s,t) = dsi ds2-, 2T7 + •- (19) 



47r2 Jo Jq {si - Pi){s2 - P2Y 

where the neglected pieces involve a complex part of the contour (the background contribu- 
tion). 

It has been shown in ref. |@] that the leading perturbative spectral function can be ob- 
tained, for this diagram, with the conditions on s and t given by (p!?]), by the 3-cut integral 
(see Fig. lb) 

.(p^ri..v^^o)^/^. ^^'''^)^^f (20) 

J [pi-k + qiY 

The evaluation of (|3) has also been discussed in p| and the answer, for the associated spectral 
density, turns out to be rather simple 



A(si,S2,s,t) = {-2Tiifj{pl,pl,s,t) 



St 



(21) 



Compared to the form factor case (see eq. ([T3|)), the discontinuity along the cut, A(si, S2, s, t), 
as given by (PTF) is Si and S2 independent. This results seems to be oversimplified. However 



such a simplification is due to the fact that 1) we have neglected the quark masses in eq. (|16[) 
and 2) we have not considered other subleading cuts in the in the evaluation of ( pTD which 
are suppressed by power of s,t or u, compared to the leading result (the single dispersive 
contribution) . 

As we are going to show next, even in the simpler conditions of a kinematics of Compton 
type, a dispersion relation in si and S2 has unobvious features, not evident from eq. (|2l|) . 

In fact, while for a vertex function there are only 3 kind of cuts and 3 independent variables 
(say Si, S2 and t), one can vary any of them and keep the others below their thresholds, 
obtaining a dispersion integral extended up to infinity, in the case of the box diagram we have 
7 cuts (respectively in s, t, u, plus the 4 external virtualities) and only 6 independent variables. 
The maximum number of variables which can be kept fixed is five. 
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A variation in the 6-th variable will then automatically affect the 7-th one. In our case, in 
particular, only 5 variables play a role since 2 external lines are fixed to be massless (photons 
of q\ = ^2 ~ 0) from the very beginning. 

Let's suppose that we fix s, t and S2 and vary si. Then the 5-th variable u = si + S2 — s — t 
can cross the threshold {u = 0) (see Fig. 2), for a sufficiently large si, and become positive. 

The diagram corresponding to this threshold is shown in Fig. 3. For m = and at fixed 
angle, however, such discontinuity doesn't play any role, as far as t is negative. 

The proof of this last statement goes as follows. In Fig. 3 (for vanishing m) the only 
possibility of having lines 1, 2 and 3 on shell at the same time, giving a non zero contribution 
to this diagram, is when these lines are collinear. For the same reason, lines 3, 4 and 5 must be 
collinear too. Therefore 1 and 5 must be collinear (i.e. qi-q2 = 0), which implies that we must 
restrict ourselves to consider forward scattering from the very beginning. The requirement of 
working at fixed (and non forward or backward) angle {—t/s fixed) allows us to exclude such 
configurations. 

A rapid look at the leading spectral density given by eq. (|2T|) shows that this is indeed 
consistent with the expression that we obtained before for A(si, S2, s, t), since at t = our 
spectral density becomes singular. 

Therefore, as far as t remains in the physical region (and m = 0), the contour representation 
of the box diagram eq. (|TBp (Fig. 2) can be sent to infinity and the Cauchy integral replaced 
by its double discontinuity eq. (|2T|), integrated over the positive semiaxis of Si and S2 

I4 = --r^ / / dsids2-, oTT 2V + subtr. (22) 

47r2 7o Jo {si - pi){s2 - P2) 

This unexpectedly simple result is a consequence of the kinematics chosen in the derivation 
of the sum rule. This result, as we are going to see, is no longer valid when an explicit mass 
dependence in the diagrams is introduced. 

Therefore, in the case m = and with on shell photons, it is also easy to write down the 
explicit expression of the dispersion integral in eq. (521) as simply as 



T{pI, pI, s, t, /i) = — — ^ -Log — ^ ) + subtr., (23) 

where we have chosen the same subtraction point fi for both momenta pi and p2 and neglected 
single dispersive parts in the integral. It is important to observe that the validity of a dispersion 
relation of the form given by eq. ( ^2] ) is crucial in order to proceed to the evaluation of the 
spectral density by Borel techniques. 

In order to illustrate how a bound on the validity of the dispersion relation comes into 
play, we have to reconsider the lowest order spectral density in its more general form, this 
time with a non zero mass in it. Below we are going to prove 3 basic facts: 

1) for a non zero mass m the spectral density developes special singularities in the complex 
Si and S2 planes; 
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2) the manifold spanned by the momenta which make the spectral density singular is a 
Landau surface for a -u-channel singularity; 

3) in the limit of vanishing mass the surface is degenerate and the spectral density remains 
regular, as shown above in eq. (pij). 

Let's proceed from point 1) and let's first discuss how a mass term is going to bound the 
region of analyticity. 

The leading spectral function - for a non vanishing m - is obtained, in this more general 
case, from the 3-particle cut integral 

Aplpl,s,t,m) = f d^k^-^^' - "^'^^-^^P' - f - - - "^'^ , (24) 



{pi — k + qiY — m? 
and it can be expressed in the form 

A(si,S2, s,t,m) = (-27ri)^J(si,S2,s,t,m) 



47rS , , 

(25) 



—t['irn?SiS2 — ^m?ss2 — Arn^ssi + Am?s'^ + Am?st — sH) 
which reproduces eq. (^Tl) in the m = case. 



The evaluation of (|2^) is not obvious and can be performed following the steps discussed 
in Appendix A, where we discuss a similar integral. For two different masses, say mi and 
m2, in the box diagram, an expression of A(si, S2, s, t, mi, 177.2) in terms of s and t, as nice as 
eq. (p5D, is very difficult to obtain. In general, rational powers of the Mandelstam function 
A(si, S2, t) will appear and cannot be eliminated by a simple reshuffiing of the indipendendent 
variables si, S2, s, and t. 

Differently from eq. (PT| ) which is mass independent and globally a regular function of 
the two pion' virtualities, eq. (p5|) shows, from its denominator, that the spectral function 
developes a singularity at S2 dependent positions in the si plane (and viceversa), as it is 
expected in general. 

We intend to elaborate some more on the origin of this singularity which is of w-type, since 
it is generated when all the 4 internal lines in the massive box diagram go on shell (see the 
discussion in 0). 

We remind here that the conditions qf, ^2 > has to be satisfied, in general, in order to get 
contribution from such discontinuity (for the box diagram), since a massless photon cannot 
decay into two massless particles at t < 0. The fact that such singularity reappears, also in 
the approximation of prompt photon emission {qf = = 0) imposed from the very first stage 
of our calculations, shows that dispersion relations in two variables, for 4-point functions, are 
very difficult to handle on a very general base even for the simple box diagram. 

Neverthless we can easily prove that A(si, S2, s, t, m), as given by eq. (^5]) , becomes singular 
when the invariants Si and S2 satisfy the Landau equation for obtaining a singularity of w-type. 

The proof goes as follows. The leading singularity of the box diagram can be described ge- 
ometrically by a simple dual diagram (generated by the momenta of the original box diagram) 



11 



|Tl|| whose volume is imposed to be zero. This condition can be usually expressed in a more 
simple form in terms of some canonical variables i/ij, proportional to the scalar product of the 
internal momenta in the diagram. In our case (for equal internal masses m) such variables are 



t - Ani^ 

2/12 = -1 2/13 = 

Si — 2im? 

2/23 = -1 2/14 = 

s — 2m? S2 — 2m? 

y'' = ^2^ y'' = ^2^- ^'^^ 

The equation of the Landau surface corresponding to Fig. 3 is in general given by a 



polynomial in yij |Tl 



1 - 2/12^ - 2/13^ - Via - 22/i22/i32/23 

-2/23^ + 2/14^2/23^ - 2?/i22/l42/24 - 2?/i3?/i42/232/24 " 2/24^ + 2/13^2/24^ 

-22/132/142/34 - 2^122/142/232/34 " 2?/i22/132/242/34 - 22/232/242/34 " 2/34^ + 2/12^2/34^ = (27) 

Inserting (^) into ( P?] ) one gets the expression 

^—{Am^SiS2 - 4m^ss2 - 4m^ssi + Am^s^ + Am^st - sH) = (28) 



which clearly coincides with the singularity of the massive spectral density in eq. (pSf ). 

It is important to realize, however, that such singularity is far enough from the region of 
analyticity in which the dispersion relations are enforced. 

In fact, the leading behaviour of such densities, even for non vanishing m, is still given 
(in eq. ^) by —l/{2st), at large s and t, and at fixed value of their ratio —t/s. It's easy to 
show from the same equation that such a threshold is located at the typical values Si,S2 ~ 
s(l + 0{m^ / s)) of the two virtualities, reasonably far away from the border of the contour of 
Fig. 2, for a dispersion relation in si and S2 to be valid. 

Whence we can summarize our findings in this section as follows 

1) We have explicitely described (by eq. (pS])) the location of the u channel cut in the 
spectral density in the two complex planes of si and S2 and 2) we have shown that the only 
parameter which controls the position of this singularity is the mass in the correlation function. 

At large s and t the term —s^ t in eq. (pSj) dominates compared to all the other terms, 
even for values of the virtualities large compared to the mass m in the correlator, but smaller 
compared to s and t. The limit m — can be studied from eq. (|2^) by sending to zero the 
mass in the numerator of such equation independently from the rest. The singularity surfaces 
degenerates, in this limit, the spectral density becomes regular and we reobtain eq. (pi]). In 



the next section we are going to extend our reasoning to the diagrams of the power corrections. 
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4 Power Corrections 



In order to analize the analyticity properties of the diagrams which appear in power corrections 
of gluonic type (see Fig. 4) we first review our method of calculation of such terms and show 
that in the massless case, again, the u = threshold, present in these diagrams, has a location 
in the complex planes of si and S2 which is parametrically controlled by the mass m of the 
OPE and disappears in the massless limit. We are then able to write down a dispersive 
representation of the diagrams for the power corrections which can be extended for si and 
S2 positive, ranging up to infinity. The pattern is similar to what already discussed in the 
previous section. Let's first observe that the power corrections can be obtained by the insertion 
of a modified quark propagator in momentum space [|l^ of the form 

I 7r^(G'^)m p{m+ ^) ^^^^ 

which has quartic powers of momenta at the denominators. (G^) is the gluonic condensate. 
Insertion of this propagator results in diagrams of the form shown in Fig. 4, for which we need 
an explicit spectral representation 

J(plpls,t)= [ 5 — 5 (30) 

dsi [ ds2-, ^'^^2w^' ^ 2\ - (2-'-) 



47r2 7o Jo (si -Pi)(s2 -pi)' 

We are interested in retrieving the explicit expression of A(si, S2, s,t) for these diagrams. 

For this purpose, as discussed in it is convenient to introduce two auxiliary masses mi 
and m2 and consider the more general spectral representation of the following amplitude 



J{pi,pl,s,t,mi,m2) 

d^k 



(pi — ky(k^ — mf)2 ((p2 — k)"^ — m^)^ (pi — k + gi)^ 

^ J{pl,phs,t,mi,m2), (32) 



dm\ dm2 

where we have defined (see Fig. 4 a) 

J{pIpI s, t, mi, m2) = / (^2 _ ^2) ((^^ _ ^)2 _ ^2) _ j^y^^^ - k + q.f ^^^^ 
and factorized two derivatives respect to the masses mi and m2 in eq. 
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The basic idea of the method used in p for the calculation of such corrections is to first 
approximate eq. ( |33D by its double dispersive part 

J{Pi,P2^s,t,mi,m2) = dsi ds2— 2T7 2^ (34) 

A-n^ Jo Jo [Si - Pi){s2 - P2) 

with a spectral density given by 

= i-2.^r I d% ^^ik^--^l)^^iiP^-kr-f2)SAiPi-m_ (35) 
J {pi-k + qf) 

is the radius of a finite contour in each plane of si and S2- Then, the dispersive part of the 
scalar integral associated to eq. ( pTD is obtained by setting mi and m2 to zero according to 
the relation 



^{si,S2,s,t) = { g^2g^2 ^''°'"'(^i>^2,g,t,mi,m2)) . (36) 

12 / mi, 7712=0 



Evaluation of (^) gives 



with 



A(si, S2, s, t, mi, ms) = (-27ri)^— — r (37) 

2dW [mi, 7X12) 



W^(mi,m2) = -Am\{Q'^f + 2mlQ'^ s + 4:{Q'^f s + 2m\Q'^ si - 2mlQ'^ si 

—2Q ssi + 2'm^Q S2 — 2Q ss2 — m^siS2 + SS1S2 (3^ 



Q"^ is given by 



= ^(Sl + S2 - t + \'^\si, S2, t)) = ^{s + u + X'/\si, S2, t)) , (39) 



with 



A(Si,S2,t) = Si + S2 + — 2sit — 2S2t — 2SiS2 

[ 2Q^ ■ 



(40) 



As described in Apppendix A, Q is the large light-cone component in the "plus" direction of 
the incoming momentum pi. Specifically, we expand all the momenta in the Sudakov base 

Pi = Qn+ + ^n", p2 = -^n+ + Qn-, (41) 
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where n"*" and n are light-cone momenta, such that n"*"^ = n ^ = 0, n+ ■ n = 0. In our 
conventions 

n+ = -^(1,0^,1)), n- = -^(l,Ox,-l). (42) 

The momentum transfer t can be expressed in terms of Q^, which can be viewed approxi- 
mately as a large "parameter" in the scattering process, by the relation 



2Q 

Notice that t = —2Q'^ for si = S2 = 0. 

Now, simple manipulation of eqs. ( P7[ ) and eq. ( PBD give 

- .3(2g2_,^)3(_2g2 + ,^)3 • (44) 

We now want to demonstrate that 

1) the only singularities which bound the dispersion integrals ( |5B| ) for the power corrections 
are of w-type 

2) they disappear in the massless limit. 

For this purpose, let's build on the results of the previous section by showing that, even 
for diagrams with gluonic insertions (Fig. 4) the singularity surface of w-type coincides with 
the singularity surface of the spectral density eq. (0). In particular, all the considerations 
made in the previous section for the case of the leading order spectral density - regarding the 
size of the region of analyticity of the correlator and the position of the w-cut - remain true 
also in this more general case. It is also important to observe that this is consistent with the 
structure of eq. ( P^ since the two derivatives which have been factorized in that equation 
(for power corrections we deal with quartic propagators) do not change the location and the 
nature of the singularities of the integral given by eq. (|33D . 

The discussion that follows also confirms the validity of the method developed in ref. 
for the evaluation of the power corrections in Compton processes. 

Specifically we are going to show that 

1) the singularity of u type in the power correction diagram (Fig. 4 b) is equivalent to 
the vanishing of W{mi, 1112) which is given in eq. (|5B[), and appears in the denominator of the 
spectral density eq. (0) of the power corrections; 

2) in the massless limit such singularity surface disappears and the spectral density (P?]) 
turns into the regular function (^4]). 

Similarly to what already discussed in the previous section, let's define the parameters Uij 
for this new diagram, as we have done in the previous section, which are helpful to describe 
such a singularity surface 

(m| + m|) t — m\ — m\ 

yi2 = — yi3 ~ 



2m^m4 2m3m2 
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The equation defining the Landau surface for diagrams of box typ, even in the presence of 
quartic propagators, is still given by eq. (|27|) . It is possible to rearrange the equation of such 
a surface in the form 



P{Q^,mi,m2,m3,m4) _ ^ 
QAm\m2rn'lm\Q'^ 

where is the mass of the i — th internal line, and P{Q^,mi) is a complicated polynomial 
in the masses mj and in the momentum Q, whose expression is too lengthy to be given here. 
Although the calculations are quite involved, neverthless it is possible to show that in the 
limit m3,m4 eq. ( ^61) simplifies since 



P{Q^, mi, m2, ms, 1714) W{mi, (47) 

and becomes 

iy(mi,m2)2 = 0, (48) 



with W{mi^ 7x12) given by eq. (^) and clearly coincides with the singularity surface of eq. 

Therefore we have shown that in the massive case - even for different masses - and, in 
particular, in the diagrams of the gluonic power correction, the spectral density developes 
singularities of w-type as conjectured in ref. 0. Such singularities are absent for vanishing rrii 
and the spectral integral, in this specific case, can be extended to the entire positive axis of 
Si and S2- Therefore, we can extend the results of ref. and write down a new dispersive 
representation of the power correction integral eq. (^) for si, S2 > as 



00 



J{pi,Pl,s,t) = --^ I dsi I ds2 , ^^^i'-^2,g,^) ^^g^j 
4vr^ Jo h (si -Pi)(s2 -P2) 

where we have omitted, again, single dispersive contributions. A check on the consistency of 
our results can be obtained immediately since A(si, S2, s,t) can be rewritten in the form 

- 327r''2(si - s)(ssi - 2siS2 - Sit - S2t + - 2i(:A^/^(si, S2, s,t) 

A(si, S2, s, t) = -— — — —J- , (50) 

[Si + S2 ~ t + A(Si, S2, s, tjjsH-^ 

and the only singularity displayed by (^) is a pole at t = 0. For scattering at fixed angle this 
pole doesn't play any role. 
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5 Radiative Corrections 



In this section we develop methods for the calculation of the radiative corrections to the lowest 
order sum rules. Although we are not going to discuss the complete evaluation of these terms 
to the actual case of pion photoproduction, which we leave as a future investigation, here we 
will describe methods which may have application to a large class of reactions of Compton 
type and are easily generalizable to the proton case. The method used is simply based on 
the observation that the complex contributions to the spectral densities can be " decomposed" 
in specific subdiagrams which involve 2- and 3-particle cuts times the lowest order result. 
After this observation, we are then able to treat these contributions in a systematic way. 
Unfortunately, from a practical viewpoint, general schemes to treat such corrections are not 
available, since double dispersion relations involve complex spectral densities beyond leading 
order. So, one has to resort on a case by case study. In particular, the treatment of the infrared 
divergences in the dispersion formalism, for simple dispersive integrals, and the problems 
which naturally come from its application, are well known |2^. Infrared divergences are often 
regulated by introducing a fictitious mass in the discontinuity function, as done, in some cases, 
for the form factor. In our case such a problem is even more complicated to deal with, since 
we are handling double dispersion integrals. As we have learned from the previous sections, 
a mass term bounds the spectral density and modifies its singularities. Particular care is 
therefore required when trying to take the massless limit in infrared sensitive integrals . 



The method that we are going to illustrate in this section has two main features, specifically 

1) it combines the good quality of the "Breit frame" (namely the fact that in the "Q"- 
variable the spectral density is polynomial) with 

2) dimensional regularization. 

These two properties are crucial in order to do the calculations in close form. Dimensional 
regularization, in particular, is used in a plane which is transversal with respect to the two 
(longitudinal) light-cone variables in the directions. The difference of the method, com- 
pared to the lowest order calculation, appears only in the angular variables, since the angular 
integral has now dimension n — 3. Given the fact that the lowest order result spectral density 
is infrared safe, the use of dimensional regularization as n goes to 4, does not modify the 
lowest order expression of the spectral density (see Appendix C). 

General arguments regarding the infrared safety of the sum rule method to Compton 
scattering have been presented in ref. 0. In a direct approach, however, one has to be able 
to show explicitely that the infrared divergences generated by the all possible cuts cancel. 
We reserve the discussion of this important point elsewhere. The treatment of the ultraviolet 
divergences, instead, is much less difficult since can be controlled by simple subtractions. In 
particular, for Borel sum rules, such difficulty is largely absent due to the fast convergence 
of the spectral densities after Borel transforms. This discussion can be addressed in a better 
way in the actual case of Compton scattering rather than in an auxiliary scalar theory, as we 
are doing it now. In particular, the use of Ward identities in a true gauge theory simplifies 
drastically the treatment and can possibly make such infrared cancellations more transparent. 
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Notice that standard tools, such as the largest time equation, commonly used in the 
calculation of single variable cuts in diagrams with any number of external lines, (see ^) are 
of very limited help in the classification of the double spectral densities of Feynman integrals. 

In order to make our discussion as clear as possible, we have decided to treat one specific 
diagram with mass regulators and the remaining ones by dimensional regularization. The 
treatment of this first diagram (which is, however, real) is quite involved, but it may serve 
illustrate quite well our general approach, and how it is possible to derive a set of basic 
rules according to which we can deduce whether a particular cut gives a real or a complex 
contribution to the OPE. 

For this purpose let's consider Fig. 5a with the momentum (and energy) flow chosen as 
shown in this figure. 

Let's first observe that the double imaginary part of such diagram (Fig. 5b) can be split into 
the convolution of a 2-particle discontinuity in the s' energy variable, specifically s' = {pi — lY 
and in a 2-to-2 scattering amplitude ( the top of the diagram), calculated at Born level, with 
on-shell external lines (Fig. 5c). 

The 2-particle cut (bottom part of Fig. 5c) contributes to the total discontinuity through 
the sub-integral 

A5''(si,S2,pi-/) 



J {{pi - kf - m?){{p2 - ky - m?) ~ 27rDV2 \ ac - hdcosO - D^^ ) 

(51) 

where 

a = s'^ — s'{si + jji^ — m?) — (/i^ — m^)si; 

i2 // I 2 2\ I / 2 2\ 

c = s — s(s2 + /U — m) + (/i — m )S2] 
d = \^l\s',fji\m') \^'\<d,S2,s')- 
qi2 = qi- q2] 

D = {ac-hdcosOf - {a^ -h'^){c^ - d^); 

cose = s'/{{s' - si){s' - S2))(2t + s' -S1-S2- S1S2/S'); 

s' = Si- 2pp- - 2p^l+. 

(52) 

and where /i and m are regulator masses. They allow us to keep under control the diver- 
gences of these sub-cuts. 
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In order to obtain the full contribution of diagram 5b we have just to convolute its top 
part, shown in Fig. 5c, with (|51|), in the form 

°° dl* r dl- I J.,^^>>i'uMt'-rny4(^n + l?-rn^)^ 
'J [Qi + — 



CO J — oo 



The integrals over the transversal components of / and on one of the two longitudinal 
components, say /~, can be done elementarily (since qi2 ■ l± = 0) while the integrals is cut 
from above and from below by the conditions on the momentum flow. The analysis of these 
conditions is algebraically very involved and is briefly sketched in Appendix B. 

It is possible to simplify eq. ( |55D even further. 

For this purpose it is crucial to observe that, in the Breit frame, eq. (|5lD has no dependence 
on the angular variables of the transversal plane, which have been introduced by the Sudakov 
decomposition of the internal momentum / 

I = l+n+ + rn- + h. (54) 

After having performed the integration over the transversal variables, the final result for 
the diagram of Fig. 5b can be expressed in the form 



'56 



where 



;„ . (56) 

2gi2 



The explicit expressions of and /i4, are given in Appendix B, together with a brief discussion. 

The contribution from diagram 5b is, however, real - for very small values of the regulator 
masses m and /i - and therefore cannot be responsible for any interference. In order to illustrate 
this last point, let's show first that the integrand in (^) is a regular function of in the 
interval of integration. Our analysis is exact in the limit of m — > 0. An independent check of 
our reasoning will also be discussed below. Let's rewrite ([5B|) as 



^1^+ + B, 



where 

Ai = 2Qi2Q^-si)iAQ^-siS2Y, 



B, = (-2g^ + si)(-2Q^ + S2)(-16g« + 16g'*s-8gV + 4g2ss^S2 



(5^ 
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Inserting this result in the denominator of eq. (|55|) we obtain 



(gr^m + qi n - 2qtqi [21+ Un -m)= ^q2(^^qa _ g^^^y ' ^ ' 



where 



iVi = (2g2-si)(4Q^-siS2f 

X (leg^ - legS + sgS^ - 4g2ssiS2 + 

iVg = m24g2(s - 2g2)(-2g2 + Si)(-2g2 + S2)(4g" - SiS2)\2Q^s - S1S2). (60) 



In the hmit m the quadratic form (pQ]) has only complex roots, since its discriminant 
vanishes. Therefore there are no singularities of the integrand of ( |55| ) in the variable A 
small finite mass term does not modify our conclusions, since it is possible to show, even in 
this case, that the two (now real) roots still lay outside the integration region in /+. 

An alternative, more direct proof of this fact can be obtained in the following way. It is 



easy to realize that any (eventual) additional imaginary part in (|55| ) is generated when the 
residual propagator of momentum pi + / (see also eq. (|53D) goes on-shell. The discussion 
outlined in section 3 shows that this is indeed not possible in the limit of a small mass m, 
unless t = 0. Therefore the integrand is regular. 

After these considerations, we are allowed to do a simple counting of the complex factors 
of 2, for this specific diagram, as a check of our conclusions: we insert of a factor (from 
the 4 single-particle cuts) times a factor of from the two independent loop integrations to 
decide of the overall factor. The result is therefore real. 

Other imaginary parts are generated by diagrams of the form 6a and 6c, in which the 
insertion of radiative corrections of vertex-type and self-energy are considered. 

Let's now analize Fig. 5d. In this case - and we shall proceed in a similar way in the other 
similar cases - we prefer to use dimensional regularization to evaluate the sub-cuts which 
appear in the relevant Feynman diagrams. The divergences, in fact, emerge at this level. 
The contribution of these sub-diagrams are then folded with a 3-particle cut integral which is 
further evaluated by using Sudakov variables in the Breit frame of the two pion' lines. 

The method is particularly convenient in order to check the infrared safety of the final sum 
of all the diagrams. 

For this purpose let's decompose 5d as shown in 5e, where the top part now involves a 
2-particle cut with 4 on-shell external lines, s' = [pi — k + qi)'^, in this particular case, denotes 
the invariant energy in the "s" channel of this subdiagram. In n = 4 spacetime dimensions 
this 2-particle cut is infrared divergent. The situation is similar to the case discussed above 
(diagram 5b), when the two masses of this diagram were set to zero. Notice that, in principle, 
also 5b can be regulated dimensionally, although the expression of its 2-particle contribution 
is far less obvious, since it involves an imaginary part with 2 off-shell external lines. 
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Let's then define the n-dimensional expression of this 2-particle cut contribution when all 
the external hues are massless 



There exist various possible way to evaluate (pl|), the simplest probably being the one 



(61) 



discussed in [0. This same method can be simply used to regulate by dimensional regular- 
ization also eq. In |]T3| , ( pT] ) is calculated by the method of Feynman parameters combined 
with the Kummer series. A simpler way to obtain this result is, as usual, to simplify the 
constraints imposed by the Dirac's delta function in a suitable frame. The derivation of these 
discontinuities is quite straightforward and we omit any of the details and quote the final 
result. We get, for the box diagram (top of Fig. 5e) 



^box 



n/2-1 



r[f-2]^ 



4 T[n/2-l]T[n~A] 
xF(l,l,^-l,l + ^)r/-^ 



(62) 



where F{a,b,c,z) denotes, as usual, the hypergeometric function [|T8]. The 2-particle cut for 
the vertex function (Fig. 6b) is instead given by 



^vertex 



{k-piY 

7r^I>/2-^^,„/2-3 

A r[n-3] * 



(63) 



The expansion around = 4 + e of eq. (^) generates a single pole in 1/e times a combi- 
nation of polylogarithmic functions, the last ones obtained from the Taylor expansion of the 
hypergeometric function 

' 'l + \^'Lh(x) + - U,(x)) + Olf")) , (64) 
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where x = 1 + t/ s' . 

The definitions of Li2{x) and 5*1^2 are given in ref. pO[. We can easily convolute this result 
with the expression of the 3-particle cut given in eq. (|20| ) to obtain 



A5, = j rk64k')64{pi-ky)64{p2-ky) A,ox{{pi - k + qf)) 
= V,j d^k5^{e)5+{{pi-kf)5+{{p2-kf) 

\2\ (/ 1 , x2\"/2-4 



X 



F (l, 1, n/2 - 1, 1 + t/{pi -k + qif) ((pi - k + qi] 
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where we have set 

r = vr"/^-^ T[n/2-2f 
'~ 2 V[n/2-l]V[n-A\ ^ ' 



Notice that all the infrared sensitive parts in ([65| ) are only contained in the factor Fi 



This is a nice feature of Dimensional Regularization, which allows us to evaluate completely 
eq. (|65D till the last stage. In fact, defining e = n — 4, the pole contribution and the finite 
parts of (pH) can be obtained from the relation 



TT 



where we have introduced a renormalization mass scale /x. Notice that it is necessary to 
perform also the remaining k integral in 4 + e dimensions, since it contributes to the finite 
part. The basic integral which appears here is 

a„=/d»M,(*')MPl^^)!]Ml|zM, (67) 

J {pi-k + gi)2 

(68) 

whose evaluation is similar to the one discussed in Appendix C. Expanding in e the logarithmic 
part of (|6^) we get 

A"' = f,T,ao{l-e/2Log[^]). (69) 
Inserting the explicit expression of ctq we finally obtain 

^ e^n-,^n-,/2 1^2 - 2]^r[n/2 - 3/2] _ 
^ r[n/2- l]r[n-4]r[n-3] ^ ' ' ' ' ' ' > 

x^<«^(l-iL0Sl=^l). (70) 



We have defined 



A = s - si- 2k'^ql - 2k ql 

(2Q^ - si)(2g^ - S2){AQ\s - + SS1S2] 

(4Q4 _ s^S2) 

B = 2\k±\\q±\ 

{2QhiS2{s - 2Q^){2Qh - SiS2]^^^^ 

{2Q^ - si){2Q^ - S2) 
(4g4 _ s,s2) ' 



(71) 



(72) 



22 



and set v = B/A. Notice that A > 0, and B > 0). Notice also that the momenta 

are fixed by the (lowest order) 3-particle cut (Fig. lb) in terms of Q, s and of the two 
virtualities Si and S2, and that f ~ l/Q^ << 1, so an expansion up to order is a very 
good approximation to the integrals above. 

The evaluation of 6a proceeds in a similar way, though it is slightly more complex. 

In this case we get 

J {{Pi-k + qiy) 

Using the explicit expression of eq. (|6^) we get 



J {{Pi-k + qiff \ 2 I, 



-k + qi] 



We have defined 



From the k integral we obtain two main terms: 



and the logarithmic contribution 



^ _ / , Mk^)s4{p, - kr)5^{{p, - kf) f {p,-k + q,y , 



(74) 



"^"/^^ {{p.-k + q,yy ^^'^ 



Their evaluation is discussed in Appendix C and their contribution is finite. The single pole 
in 1/e, in diagrams of this type, is generated only by the 2-particle cut of the vertex function 
shown in Fig. 6b 



ur. „ „ ^,n-5/2 r[n/2-2] 



A^" = ^J(g,Si,S2)7r 



^' r[n-3]r[n/2-3/2] 

. (t-^'-^I^ F(1.3/2. „/2- V„)(1 - Io,|(l±^M,) 

\ 1 — 3J //^ 



V 



-TX e 



Again, in practical applications, just the first few terms in v are necessary. 



(75 
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A radiative correction of another type, also responsible of the generation of imaginary 
parts is shown in Fig. (6c), and it corresponds to a self energy insertion on the fermionic line 
at the top of the box diagram. An odd number of particle cuts is involved in this diagram. 
As we have discussed in the former cases, also here, again, we can unfold the 5-particle cut 
into a 2- times a 3-particle cut. Both cuts are regular in n = 4 dimension, so there is no need 
to use dimensional regularization. For the 2-particle cut we get {q is the external momentum) 



It is not so difficult to show that the condition {pi — k + gi)^ > is identically satisfied for 
all possible values of k in the integral. 

Evaluating this integral at n = 4 we get 

2stA{l - v"") ^ 

The evaluations presented in this section, as we have seen, can be carried on up to the 
last stage. This is a remakable feature. A simple factorized structure for the imaginary 
parts generated by the radiative corrections to the lowest order sum rules emerges. All these 
corrections, (see Fig. 7), can be easily evaluated by these methods. 

It is not difficult to realize that the set displayed in Fig. 7 is the complete one. The 
diagrams shown in Fig. 8, for instance, have not been added to the fist. Let's see why. 

As an example, let's consider the diagram shown in Fig. 8a. From a first inspection, it 
seems that this diagram should be included to the list since the number of its cuts is odd. 
However, it is easy to realize that it involves a subcut of t-type, which is non vanishing only 
if t is in the unphysical region {t > 0). Therefore 8a can not contribute as far as t remains 
in the physics region. 

Other contribution to the real part of the spectral density are dispalyed in Figs. 8b and 
8c. Each of these two diagrams has a maximal number of internal on-shell lines and no other 
additional line can go on shell. In fact any other additional line, if on-shell, induces at least 
one vertex with 3 on-shell massless lines in the new diagram. The presence of any vertex of 
this type then, will force all the 3 particles at that vertex to be coUinear. This condition is, 
in general, too strong to give a non vanishing discontinuity, except at specific points in the 
phase space of the process. The discussion presented in section 3 may serve well illustrate this 
aspects. In fact, in that case, t was forced to be zero in the presence of such configurations in 
the reduced diagram. 




(79) 



and the 5-particle cut can then be set in the form 




(80) 
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6 Conclusions 



We have discussed in detail some crucial aspects concerning the analyticity properties of the 
spectral densities in Compton processes. The analytical bounds on the region in which a 
dispersion relation is valid is controlled by the quark mass in the expansion of the correlator. 

For nonvanishing m the spectral density has special singularities, absent in the form factor 
case, which have been identified as Landau surface of w-type. At order as they are important 
in the analysis of the radiative corrections. 

As we have seen, in the massless simpified picture of such corrections emerges. 

Remarkably, the evaluation of the complex parts, responsible of effect of interference in the 
sum rules, can be done in close form. For this purpose, Sudakov methods in n dimensions 
and dimensional regularizations have been developed and combined in an original way. Since 
it is generally believed that, at intermediate energy, the phases of elastic scatterings are of 
non-perturbative origins, it is interesting to see how relevant the "Feynman mechanism" 
which motivates the sum rule description is - compared to hard scattering factorization - in 
the description of such effects. 

We reserve to discuss the phenomenological application of our results to the evaluation of 
the phases of pion photoproduction elsewhere [jlTl . 

However, from this work, a well defined structure for the complex spectral densities for 7 
7 collisions emerges, which can be enforced in an actual calculation. 
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Appendix A. The light-cone formulation. 

Here we summarize briefly the kinematics of Compton processes and point out the impor- 
tance to work out the expression for the spectral densities in the light-cone frame since, in 
such a frame, they become polynomial. The simplification is enormous since the Mandelstam 
function (see eq. 0) -which appears in the form factor case- is absorbed into the definition of a 
large "plus" ligh-cone momentum of one of the off shell parton' lines. It is important to notice 
that, since from the very beginning we have been working with a diagrammatic expansion of 
current correlators, in each diagram in the O.P.E. (see Fig. 1) the virtuality of the pion' lines 
will, in general, be different. 
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Nevertheless we can picture Fig. 1 as an ordinary scattering process and select a special 
frame in which all the transversal momentum belongs to just a single particle (the photons). 
Let's define 



s = {pi + qif t={q2-qi)\u = {p2-qi)\ (A.l) 

with 

s + t + u = si + S2. (A.2) 
We also expand the momenta of the incoming photon as 

qi = + gfn" + gi_L. (A. 3) 

In the Breit frame of the incoming meson we have 

u = iP2- qif = - ^ + (A.4) 
Covariant expressions for qf and can be easily obtained in the form 

In the light cone frame the spectral function can be calculated to lowest order in terms of 
3-cut integrals of the form [0 



[pi - k + qiY 



I'[f{k\ k ■ Pi, ...)] = j d'kf{k\ k.p,, ...)6+ie)64ip^ - k)')64ip2 - kf), (A.8) 

and as in standard evaluation of Feynman diagrams, we can reduce tensor discontinuities to 
scalar one by any reduction procedure. 

Mass dependent integrals can be analized quite easily in the light cone frame. 

For illustrative purposes consider 



</34(Sl,S2,S,t,mi,m2 

{pi-k + qiY 



^,, Mk')S+{{P2 - ky - ml)5^{{p, - kf - ml) _ ^^ g^ 
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In this integral, the components of the internal momentum are fixed at the values 



- siS2 



(A.IO) 
(A.ll) 



We get 



with 



ang 



2^ rf^ 



A + Bcos[e] 

TT 



2 -r (A.13) 



A2 = s + mi-2k'^p^ -2k p^-2k^q^ -2k q^; 

B2 = (Sqtqi{2k+k- -mi)f''^ . (A.14) 



Remarkably Tang is a rational function of S2, t) in the Breit frame even for nonvanishing 

mi and m2. This is due to the identity 

{-AQ^ + s,S2f 

w(mi,Q^si,S2) = (-AniiQ^ + 2m2Q'^s + AQ'^s + 2miQ'^si 

—2m2Q'^si — 2Q'^ssi — 2miQ^S2 — 2Q'^ss2 

-miSiS2 + SS1S2 )^ (A.16) 



which gives 



TT 

^34(si,S2,S,t) = — — -, (A.17) 

26W{mi,m2) 



and is the result given in eq. (pH]). 
Appendix B. Radiative Corrections 
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The evaluation of eq. ( ^31) proceeds as follows. The integrations in the / , and in the 
angular variable in the transversal plane can be done explicitely. We define 

cos (f) = 1^"*" l'^] , (B.l; 

and since qi± = q2±, we obtain the relation qi2± -1 = 0. The only angular dependence in 0, in 
the integral, is contained in the factor {qi — /)^. Therefore we obtain 



'56 



1 /■27r Aof/+ / 



dr 



^ql^Jao Jo {{qtl'm + qil^ -2\q±\\l±\cos(, 
4gi2 Jao ( (q+i^ ^ _ 4q+q-l+l^y^^ ' 



which gives, after integration, eq. (|55D. 

Ill and fi2 are momenta at the boundary, fixed by the conditions on the energy flow. In 
order to fix these two momenta let's rewrite these conditions in the light cone variables 



1+ +lm> 0, 

qt + 1^ + qi + r >0, 

qt2 + 1^ + + r > 0, 

s' = (pi - If = si- 2ptr - 2p^/+ > 0. (B.2) 

In the last condition above, which is simply the requirement of s' > 0, we have for simplicity 
neglected the contributions coming from the regulator mass /i and from m (specifically {s' > 
(m + /i)^)). This approximation is not going to affect our conclusions in any significant way. 

The analysis of these conditions is quite involved. Here we summarize the basic results. 

In the Breit frame these conditions are polynomial constraints in the variable and can 
be re-stated respectively in the form 

/+ > /Ui, 
/+ > /U2, 
> /^3, 

(B.3) 

where 

^ {-2Q^ + Si){-2Q^ + S2) 

2Q(4g4 _ siS2) (-4g4 + siS2)2(16g8 - leges + sQ^s^ - 4Q^ssiS2 + slsl) ' 

(B.4) 
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= 

Qjs - 2Q^){-2Q^ + S2)(-16g6 + SQS + 4QSi - AQhsi + 4gS2 + sls2 - sisj) 

(-4Q2 + s^ + S2)(-4Q4 + siS2)2 



(B.5) 



^ (-2g^ + S2)P3(Q,gl,g2,g,t) , ^ 



and where 



A(<5,Si,S2,s,t) = (-32g^° + 32Q^s - leg^^s^ - 16QSsi + 8QVsi 

+16g^S2 + 8Q^ssiS2 + 8g^s?S2 - 4g2ss2^2 
-8Q^sisl-2Q^slsl + slsl), 



(B.7) 



= 2g^(-2g2 + s2)(i6g^-i6g^s + 8g^s^ + 8g^ssi - 

AQ'^s^si - AQ^ssiS2 - 4:Q^sls2 + 2ssls2 + sjsl). (B.8) 
It is possible to prove the vahdity of the following inequalities 

< IJ,2 < fi3 < fii < < fii (B.9) 

and hence the integration over the variable in (p.2| ) is cut from below by fii and from above 
by /i4. 



Appendix C. Sudakov methods in n spacetime dimensions 

In this appendix we cover some of the technical derivation of the n dimensional disconti- 
nuities, obtained by using a special technique which is essentially based on the use of n — 2 
transversal dimensions for the Sudakov variables in the Breit frame [|l^ . 

The evaluation of the longitudinal parts of the integrals proceeds in the same way as 
illustrated in Appendix B. Therefore the components of the momenta inside each loop of 
integration are fixed. The evaluation of these longitudinal parts brings in a Jacobean, which 
is polynomial in the momentum Q, which is evaluated on the two dimensional subspace of 
n^. The transversal integral, instead, is expressed in terms of n — 2 angular variables. Of 
these, n — 3 can be integrated trivially, while the remaining one, which we call 6, is non trivial. 
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The angular integral so generated are finite, but we need there expression up to 0(e), since 
also these terms contribute to the finite part of the spectral densities. We illustrate here the 
derivation of ai and (J2, as defined by eqs. ( [76D and ([77|) . 
We obtain ^ 

CTi = —J{Q,Si,S2)ujin)Tang,iin,v) (C.l) 

where 

JiQ,SuS,) = -—^ V (C.2) 

2[4:Q'^ - S1S2) 

is the two-dimensional Jacobean, while 

^n/2-3/2 

uj(n) = 2—- — (C.3) 

^ ' r[n/2-3/2] ^ ^ 

is a volume factor obtained from integration over n — 3 angular dimensions. The nontrivial 
part is contained in 



/•ZTT Sin U 

Jo (1 + vcos oy 

which is of hypergeometric form 

r„„,,(„. . ^imzl^F(2, „/2 - 3/2. n - 3. (C.5, 
Using the quadratic (Gauss) relation 

F(a, 6, 26, z) = {l- ^-)-'^F{l/2a, a/2 + 1/2, b + 1/2, (^)') (C.6) 
we can re-express ( |C.5| ) in the form 

T^ngM = 2n-5 rK2 - 3/2]^ ^ _ 
r[n — 3J 

Notice that, being u < < 1, the hypergeometric function in ( |C.7| ) can be approximated, for 
all the practical purposes, by the first 2 terms, namely up to order v. 

(72 can not be reduced to a complete explicit expression, because of the logarithmic cor- 
rections, which are not exactly of the hypergeometric form. This contribution can be directly 
evaluated with n set to be 4 from the beginning, since these corrections already appear with 
a factor of e away from n = 4. 

For this second integral, after that the longitudinal integration has been performed, we 
obtain 

J(Q^Sl,S2) / / X ON 

0"2 = Uj{n)Tang,2{v,n), (C.8) 
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where 



(1 + vcos oy 

Log[—]Tang,l{'^,v) + Tang ,3{n , v) . 



(C.9) 



We have defined 



on— 5 /•! 

Tang,3^J^-^ tLo^? [(1 + t;) (1 - ^ t)]^"/^"^^/^ ( 1 - tf/^-^/^ (1 - ^ ^)-^ (C.IO) 



[l + vyjo 

and z = 2f /(I + v) ^ 2v. Expanding the logarithmic part to lowest order as Log[l — zt] ^ 
—zt + 0{v^) and taking n = 4, we get 

2v 

Tang,3-Log[l + v] T,„,,i (4, t;) + -^^^5 [3/2, 1/2] + Oiv^) (C.ll) 

Combining all the terms together we finally obtain 

W = Lo<7[^^i±^]T,„,,i(4,t;) + -^+ 0{v'). (C.12) 
/i^ [l + v) 

Using these results for the angular integrals, it is straightforward to obtain the explici 
expressions of A^*^ and A^°. 



Appendix D. Borel methods 

In this appendix we present a discussion of how and of when the Borel method [jl| can 
be applied to Compton scattering. We have organized the material starting from 2 point 
functions. Then we come to the pion form factor A simple way to invert Borel transformed 
amplitudes for 3-point functions is presented. Then we show that the method can be extended 
to lowest order Compton scattering. In this case it is shown that the divergence of the 
Schwinger representation in the box diagram can be kept under control by working in the 
euclidean region. 

To make our discussion self-contained we first illustrate briefiy the simplest application of 
this technique to correlators of lower rank. 

Let's consider the dispersion relation of the polarization operator 



^, 2x 1 f°° A(s)ds , . ,^ ^, 

ll(gi) = — 7T- + subtractions, (-U-lj 

TT 7o (s - gf ) 
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with a singularity cut starting at = 0. Eq. (|D.1| ) can also be written in the form 

1 fOO PCX) „ 

n(g2) = _/ ds rfaA(s)e-"(^-''?\ < g, (D.2) 

TT JO JO 

where we have used the exponential parametrization of the denominator. 

The Borel transform in one variable is defined in its differential version by the operator 



Q2/n=Af2(n- 1)! dQ' 

denotes the Borel mass. 
It satisfies the identity 

B{Q^ M2)e-°^' = 5(1 - aM^) aQ^ > 0. (D.4) 

Acting on the polarization operator Il{qf) with the Borel transform we get the usual 
exponential suppression of the higher states 

M^B{-ql ^ M2)n(g2) = _ / ds A{s)e-''^' . (D.5) 

TT Jo 

At this point we can Borel transform once again eq. (p.5| ), with respect to the inverse Borel 
mass in order to obtain 

B{l/M^ u){M^B{-ql ^ M^)Yi{ql)) = ^A(l/z/). (D.6) 

Therefore, by acting iteratively with Borel transforms on n(q'i), we obtain an expression from 
which we can easily identify the spectral weight A(s) of eq. (|D.1| ). 
Let's now turn our attention to vertex functions (see also Fig. Ic). 

The Borel transformed amplitude for the pion form factor (the complete spectral density) 
is given by il]] 



1 roo roo 

mlMlq') = - dsJ ds,plT{s^,S2,q 

TT^ JO JO 



2^ „-si/M2-S2/M. 



2 



27r2(M2 + M|; 



dxxil 



X) exp 



-xq 



(1 -x)(Mf + M| 



(D.7) 



To isolate the spectral function from eq. (|D.7D we need to use Borel transforms and act on it 
with the differential operator 



B{l/M^ l/vi)B{l/Ml l/v2)MlMl. 



(D.8) 
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Here we discuss a possible way of doing this. The inversion of (p.7| ) can be obtained by taking 
inverse Lapace transforms twice - respect to - of this equation. In fact, for a given 

function F(M^), the following identity 



L-\llM' ^ u)F{M') = {l/u)B{l/M' l/u)F{M') 



relates the differential operator given in ( D.3 ) to the inverse Laplace transform 

^ rc+ioD 



L-\l/M^ -^u) = — / d{l/M^) exp (v/M^) . 

2,TTi Jc~ioo ^ 



(D.9) 



(D.IO) 



Defining = /ii, 1/M| = /i2, then , we are required to act with the operator defined in 

eq. (|D.8| ) on the integral function 

1 



x(/^i,Ai2,g 



dx 



X 



x + iy 



■ exp -xQ^/i2 + 



2 2 



Hi + /i2 Jo 
Using the gaussian relation 

/OO 
da exp{—ka'^ — aa 
-OO 

we can rewrite eq. into the form 



1^1+ 1^-2, 



(D.ii: 



(D.12) 



do' / dx — z — — z — TTT^ -. (D.13) 



By using the relation 

on (p.l3|) we finally get 

_ /-OO a^t^a; .^FT exp((-sgV2-2xi/2g/i2a-/i2)«^i 

= L (/i2,z^2) / 7 — ^T7 / da 



{x + iyJ-^ (yi-a^yl^ 
, (-cr + /i/2)V 

^ (a_/l/2)2 + ^2)4^1/2(^^ _^2)l/2 

{\l^)q\dldq^fY{q\vi,v^), (D.15) 
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where 

/ = - Z/i - Z/2. (D.16) 

By redefining Ui = Si, i = 1,2, and after doing the exphcit integration of in eq. ( p.l6| ), it is 

easy to relate this last result to the 3-particle cut integral for the triangle diagram (Fig. lb), 
in the scalar case, 

si, S2) \ f +\ m 17^ 
^ = A3(si,S2,t) (D.17) 

The derivative with respect to qf in eq. (p.l5|) takes into account the fermionic character of 
the propagators in Fig. lb compared to the scalar case (see eq. ( p.l7|) ). The spectral function 
for the pion form factor can then be expressed in the form H] 

PwT^^iSl, S2, t) 

1 

fD.18) 




{{Sl + S2-tf-^S,S2Yl^ 

with t = —q^. 

For four-point functions this procedure simplifies considerably. The method can be applied 
exactly as in the form factor case, although it is necessary to work in the euclidean region from 
the beginning. However 1) the amplitude has first to be continued in the unphysical s and t 
region before any Borel transform is taken 2) afterwards the result has to be continued back to 
the physical region. These two steps allow us to bypass problems related to the divergences of 
the Schwinger representation (see the discussion below) of the amphtude. Since, by doing so, 
we are comparing timelike and spacelike representations of the same amplitude, we analize 
step by step the effect of the Borel transforms both on the timelike amplitude and on the 
spacelieke one. As we have already emphasized, the method works only if the dispersion 
relation is unbound, i.e. the dispersion integral can be extended to infinity for both si and 
S2- We have pointed out that for massless correlators this is indeed the case. 

Let's work in the euclidean region of Ti{p\^p\, s,t) with spacelike external invariants 
{PiiVi < 0, s = (pi + qiY < 0, t = {q2 — qi)^ < 0). It is then possible to relate r4 to 
its euclidean continuation T^e and use the Schwinger parametrization for the latter: 

= J d'^k j^[dai] exp (-ak'^ - (3{p[ - kf - -f{p[ -k + q[f - e{j)'^ - kf) , (D.19) 

where ctj is a short notation for the proper time parameters. Divergences, in this representa- 
tion, reappear when we move into the physical s,t region. We perform the integration over 
the loop momentum in eq. ( p. 19 ) to get 



Z4E = / dxidx2dx'idxi5{l — xi — X2 — x^ — Xi) / T^dT^e^ (D.20) 
Jo Jo 
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where 

r = -S(Ai(xi)si -A2(x,)s2- A3(a;,,t)) (D.21) 

with suitable expression for Ai, A2, A^. In the region where = —p'1 p\ = —p'2 and, in 
general, pi ■ pj = —p\ ■ p'^, T^e is related to T4 by an analytic continuation 

TM.Pl^.t) = iT,E{p'lp%s\t') (D.22) 
In the euclidean region the form of the dispersion relation becomes 

Exactly as in the case of the polarization operator (eqs. ( p.l| ) and ( p.2|) ), we can now apply 
Borel transforms on eq. (p.23|) to get 



r,E{MlMls',t'' 



B{p[' ^ l/M^)B{p'2' ^ l/M^)nE{Pi,p'2,s',t 



Jo Jo 

I dxi dx2 dx-i 6(1 - AiJ:M^)5(1 - A2M^i:)e'^^'i:^ dT. 
Jo 



47]-2 



(D.25) 



where we have defined 



(l-a;i-a;2) 
b\ = s— 



X2 

h2 = iz£^l±M. (D.26) 

Xi 

Notice that B{p'-^ Mf) = B{—pf — > Mf), since < by assumption. 

By acting similarly on T4 as given by (0) - with B{—pf 1/Mf) - we also get 

v{MlMls,t) = -— dsJ ds2Ais^,S2,s,t)e-'^/^'?-^-/^''- 

477^ JO Jo 

= B{-pI^MI)B{-pI^MI)UpIpIs^). (D.27) 

Reapplying the Borel transform (this time with respect to 1/Ml) sequentially on both 77 
and TjE and using eq. ( p.4|) we get respectively 

B{l/Ml ^ z/i) B{l/Ml ^ V2) {MlMl 7]e{MI s', t')) 
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1 rirf^ 

Jo Xfx2 

^ As(-l/z/i,-l/z/2,s',t') (D.28) 



for the eucUdean part, and 

B{l/Ml z/i) B{l/Ml U2) {MlMl r]{Ml, s, t)) 



^ A(l/z/i,l/z/2,s,t). (D.29) 



47r^z/iZ/2 

in the physical region. As we can see from eq. (|D.28| ), the integral on the two parameters Xi 
and X2 is now trivial and can be performed straightforwardly. A simple calculation then gives 

47r^ 

Ae{-1/i^u -1/1^2, s't') = -. (D.30) 

Using the identification 1/z/i — > si,l/z/2 — S2, analytically continuing back (the euclidean 
result in (p.3CI|) with the prescriptions s' —s t' ^ —t and using eq. ( p.22|) we finally get 

Aiplpls,t) = -At^, (D.31) 
st 

valid in the physical region. This result agrees with the one obtained from the cutting rules 
eq. (pri). The Borel method allows to overcome, in this simpler case, the difficulties related 
to the evaluation of discontinuity integrals. The original discussion of this section therefore 
provides an extension of the method to Compton scattering, from the known case of the form 
factor IjlJ. As we have extensively discussed, however, the presence of w-channel singularities 
makes it inapplicable to Compton scattering in full generality. 
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Figure Captions 



Fig. 1. (a) The Compton scattering amplitude; (b) the 3-particle cut for Compton scattering. 
One internal line is off shell, (c) The triangle singularity for the form factor to lowest order. 

Fig. 2. The integration contours for the scalar amplitude in the pf plane, which includes the 
M-channel threshold (thick black line) . 

Fig. 3. The w-channcl cut of the box diagram. 

Fig. 4. (a) The diagram for the gluonic power corrections, (b) The 3-particle cut for the 

same diagram; (c) its u channel cut. 

Fig. 5. (a) Radiative correction to the box diagram in the scalar amplitude. The line with 
the encircled "1" denoting a scalar of mass (b) A cut generated by (a), 
(c) Its decomposition as a convolution of a 2-particle cut in the s' energy (bottom part) with a 
2-to-2 on-shcll scattering amplitude, (d) The cut giving the imaginary contribution generated 

by (a) and (e) its factorized form. 

Fig. 6 (a) A subdiagram of vertex type with a 2-particle cut. (b) Its factorized form, (c) The 
self-energy 2-particle cut. (c) The imaginary contribution to the spectral density obtained by 
the insertion of (b). 

Fig. 7 The complete OPE of the imaginary parts of the double spectral densities for Compton 
scattering. 

Fig. 8 (a) A diagram containing a cut of type t. (b) The discontinuity of t type extracted 
from (a), (c) A possible (real) cut which is not factorizable. (d) The same as in (b). Notice 
that 3 propagators are off-shell. 
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